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Abstract 

We describe a simple method of including dissipation in the spin wave band structure of a 
periodic ferromagnetic composite, by solving the Landau-Lifshitz equation for the magnetization 
with the Gilbert damping term. We use this approach to calculate the band structure of square 
and triangular arrays of Ni nanocylinders embedded in an Fe host. The results show that there are 
certain bands and special directions in the Brillouin zone where the spin wave lifetime is increased 
by more than an order of magnitude above its average value. Thus, it may be possible to generate 
spin waves in such composites decay especially slowly, and propagate especially large distances, for 
certain frequencies and directions in k-space. 

PACS numbers: 



1 



The existence of a periodic superlattice strongly affects many types of excitations in solids. 
For example, the electronic band structure of a conventional semiconductor or semimetal 1|, 
and the dispersion relations of electromagnetic waves {2!, elastic waves S-gI, and spin waves 



7- 



ll| are all greatly influenced by a periodic superlattice potential. In many cases, such 



potentials can give rise to new, and even complete, electronic, photonic, elastic, or magnonic 
band gaps which may have important implications for the properties of these materials. 



These excitations have, by now, been extensively studied numerical 



ly and analytically, using 
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a variety of methods, and have been probed in many experiments 

In the present paper, we consider a particular class of such excitations, namely, spin waves 
in periodic magnetic materials. Such magnetic superlattices are often called magnonic crys- 
tals. We go beyond previous work by calculating the spin wave lifetimes in such materials. 
Our most striking finding is that the "figure of merit" (FOM) of these spin waves (product 
of spin wave frequency and lifetime) is strongly dependent on the Bloch wave vector k, even 
though, in our model, the same spin waves would have a k-independent FOM in a homoge- 
neous magnetic material. This strong k-dependence suggests that magnetization in periodic 
magnetic materials may be transported most efficiently by spin waves propagating along 
special directions in k-space. Possibly this k-dependence could be tested by experiments in 
which spin waves are launched in particular directions corresponding to the largest FOMs. 
This spin wave generation could be accomplished using real magnetic fields, or (via the spin 
torque effect IGj) using spin currents. Measurements of spin wave lifetimes might be carried 
out, e. g., by neutron spin-echo techniques/citebayrakci. 

Our calculations are carried out for an array of infinitely long circular cylinders made 
of a ferromagnetic material A embedded in another infinite ferromagnetic material B. All 
the cylinders are taken to be parallel to the z axis and their intersection with the xy plane 
forms a two-dimensional periodic lattice. We consider two arrangements of such cylinders: a 
triangular and a square superlattice. An external static magnetic field Hq is applied parallel 
to the axis of the cylinders, and both ferromagnets are assumed to be magnetized parallel 

to Hq. 

The equation of motion for this periodic composite is given by the Landau-Lifshitz-Gilbert 
(LLG) equation jl8l|: 

d a f d \ 

-M(r,t) = 7/^oM(r,t) x Heff(r,t) + — (r) f M(r,t) x ^M(r,t) j . (1) 
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Here 7 is the gyromagnetic ratio, which is assumed to be the same in both ferromagnets, 
Hcff is the effective field acting on the magnetization M(r,t), r is the position vector, a is 
the Gilbert damping parameter and is the spontaneous magnetization. For this inho- 
mogenous composite Hefj can be written as 

Heff(r, t) = Hoz + h(r, t) + ■ M(r, t), (2) 

where h(r,t) is the dynamic dipolar field and A denotes the exchange constant. The 
last term on the right-hand side of eq. (2) denotes the exchange field. For the two- 
component composite we consider, the exchange constant, the spontaneous magnetiza- 
tion and the Gilbert damping parameter take the forms A{r) = As + Q{r){AA — Ab), 
Ms(r) = Ms^B + 0(r)(^s,A — Ms^b), and Q;(r) = + 0(r)(aA — where the step 

function 0(r) = 1 if r is inside ferromagnet A, and 0(r) = otherwise. 

We separate the static and time- dependent parts of the magnetization by writing 
M(r, t) = MgZ + m(r, t), where m(r,t) = m(r)e~*'^* is the time-dependent part of the 
magnetization. The time-dependent dipolar field h(r)e~*'^*, where h(r) = — V\E'(r) and ^^(r) 
is the magnetostatic potential. Since V ■ (h(r) + m(r)) = 0, the magnetostatic potential 
\I^(r) obeys V^\I^(r) - V • m(r) = 0. 

Within the linear-magnon approximation [lol . the small terms of second order in m(r) and 
h(r) are neglected in the equation of motion. This is equivalent to setting m(r) ■ z = 0|20|. 
Substituting the above equations into eqs. ([1]), we obtain 

Ms 9^ 

iilm^{v) + V ■ [QVmy(r)] ~ "^2/(^) ~ + = 0, 

iQmy{r) -V ■ [QVm^{r)] + m^{r) + inam^{r) = 0, (3) 

IlQ OX 

where = u/{\'j\noHo) and Q = 2A/ {MsIIqHq). 

Next, using the periodicity of Q, Ms and a in the xy plane, we can expand these quantities 
in Fourier series as Q(x) = Q(x, y) = YIg Q(Gr)e'*^''', with analogous expressions for Msi;^) 
and a(x). Here x and G are two-dimensional position and reciprocal lattice vectors in 
the xy plane. The vector r = (x, z), but none of the above quantities will have any z 
dependence for the composite we consider. The inverse Fourier transforms are of the form 
Q{G) = ^ / / (i^xQ(x)e~*^''', where S is the area of the unit cell; similar expressions hold 
for Ms{G) and a(G). 
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To calculate the band structure for spin waves propagating in the xy plane, we con- 
sider the two-dimensional Bloch vector, k and use Bloch's theorem to write ma;(x) = 



ikx 



G''"x,KlG)e**^-'', rriyiyi) 



e'''-^^^mj,,K(G)e^<^-'', and ^(x) = e*'^'' ^K(G)e 



G X 



After some straightforward algebra, the equations of motion reduce to 



^fi^i(G,G') 



"^x,k(G) 
"^y,K(G) 



5^M(G,G') 



"^x,k(G') 
"^^/,K(G') 



(4) 



i(G,G') 



(5) 



the 2x2 matrix 

5gg' a(G - G'; 
-«(G-G') (5gg', 

where 5gg' is the Kronecker delta and the four components of the 2x2 matrix M(G, G' 
are given by 



M(G,G').. 
M(G, G')., 

M(G, G')y, 
M(G, G')yy 



Ho 



(K + G 



'\2 



5gg' + QiG - G')(K + G) • (K + G') + 



M,(G- G') {Ky + G: 



y> 



-Sgg' - Q{G - G')(K + G) • (K + G') - 
M,(G-G') iK., + G'JiKy + G'y) 



Ho (K + G')2 

M,(G-G') + 



(K + G'Y 



Ho 



(K + G')2 



(6) 



On left-multiplying eq. (jl]) by the inverse of the matrix A, we reduce the band structure 
problem, including Gilbert damping, to that of finding the (complex) eigenvalues of A~'^M. 
A similar plane wave expansion has been previously used to calculate the magnonic band 

n n 

structure, for the case of zero damping, by several others (see, e. g., Refs. [7[ and |18l|). 

We have used this formalism to calculate band structures for both a triangular Bravais 
lattice, with basis vectors ai = ax, a2 = a (^^x + -^yj , and a square Bravais lattice, with 
ai = ax, a2 = ay, where a is the edge of the magnonic crystal unit cell. Since Fourier 
transforms are available analytically for cylinders of circular cross section, the band structure 
is easily calculated in this plane wave representation. 

In order to solve Eq. (j4j), we restrict the sum over G' to the first 625 reciprocal lattice 
vectors, which requires the diagonalization of a 1250x1250 complex matrix. The resulting 
eigenvalues of the matrix -B(G, G') are all complex. For a given k, the imaginary part of the 
eigenvalue for gives the spin wave frequency, while the real part represents the inverse spin 
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wave lifetime. We have found that both the frequencies and hfetimes are well converged to 
within 0.1 % for this number of plane waves. 

For each eigenvalue, the figure of merit (FOM) mentioned above is the ratio of the 
imaginary part to the real part of the eigenvalue. If the Gilbert damping parameters a a = 
the FOM would be same for all k's and all bands. By contrast, when 7^ we find 
that the FOM varies from band to band and depends strongly on k. In particular, the FOM 
is particularly large in certain high symmetry directions. As a result, spin waves will have 
a longer lifetime when they are launched at special k values and with special frequencies. 

We first consider the case of zero damping. In the left panel of Fig. 1, we plot the band 
structure of a composite of Fe cylinders arranged on a triangular lattice and embedded in 
a Ni host, as calculated at an applied field ^qHq = O.IT. The lattice constant a = 10 nm 
and the Fe filling fraction / = 0.5 (f is the area fraction occupied by the cylinders). The 
center-hand panel shows a similar composite, but for Fe cylinders arranged on a square 
lattice, again with / = 0.5. The right-hand panel shows the Brillouin zones of the square 
and triangular lattices with symmetry points indicated. In calculating the band structure, 
we use an exchange constant and spontaneous magnetization at room temperature of 8.3 
pj/m and 1.71092 x 10*^ A/m for Fe, and 3.4 pj/m and 0.485423 x 10^ A/m for Ni|21|. 
We have not found band structures for exactly these materials in the literature, but when 
we carry out analogous calculations for Co cylinders in a Permalloy matrix (not shown), 
usin g th e plane wave method, we obtain nearly identical results to those found by Vasseur 
et a/[l8|, who also used a plane wave expansion. 

In Fig. 2, we show analogous calculations including damping for a square lattice. We use 
the same parameters, magnetic field, and value of / as in Fig. 1, except that the Gilbert 



damping parameters are ape = 0.0019 and aArj=0.064, following Ref. [22|. In the left panel, 
the width of each cross-hatched region is proportional to the figure of merit (FOM) for 
the given band and k value. The right panel shows the FOM for the fourth lowest spin 
wave band, as a function of magnonic crystal wave vector k, along specified directions in 
the superlattice (or magnonic crystal) Brillouin zone (SBZ), and at three different filling 
fractions /. The inset again shows the SBZ and symmetry points. We plot the first nine 
bands. The scales for the FOM and the real frequencies are different, as indicated. 

In Fig. 3, we show the corresponding quantities for a triangular magnonic crystal, again 
using a superlattice constant 10 nm and / = 0.5. The other parameters are the same as 
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in Fig. 2, except that now the right hand panel shows the FOM for the third lowest spin 
wave band. In Fig. 4, we show how the FOM for the optimal special symmetry points of 
Figs. 2 and 3 and bands depends on the superlattice filling fraction /. Note, in particular, 
that the FOM increases strongly near the close-packing values of / for both the square and 
triangular lattices. 

The most striking feature of these plots is the strong dependence of the FOM on both k 
and band index. For example, in the square superlattice, the FOM is largest in the fourth 
band at the symmetry point M, and in the triangular superlattice, it is largest for the 
third band at K. The physics behind these strong maxima in the FOM is that, in both 
superlattices, the spin waves at these k-points propagate mainly through the Fe host, which 
is the low-damping component. This result suggests some possible ways to increase the FOM 
even further at these points: if we can arrange that a spin wave propagates entirely through 
the low-dissipation material, this should give an FOM close to the theoretical maximum, 
which is that of this material in its homogeneous form. Thus, a judicious exploration of 
different periodic composites made of Fe and Ni, or other materials, could well lead to an 
even stronger dependence of spin lifetime on k value. 

We should add a few words of caution regarding the "spin waveguiding effect." In prin- 
ciple, a measure of distance traveled by apropagating spin waves is given by the coherence 
length (or spin wave mean free path)Zc 23|]. This coherence length, for a given band n at 
wave vector k, is defined as /c(k, n) = |V^„|/7kn, where V^^ represents the group velocity and 
7kn represents the imaginary part of the eigenfrequency, i. e., the inverse lifetime. Since the 
group velocity may itself depend strongly on n and k, the behavior of /c(k, n) may be quite 
different from that of the lifetime. Nevertheless, we expect that /^(k, n), like r(k, n) and the 
FOM 7kn, will depend strongly on both k and n, with sharp extrema near special symmetry 
points. Hence the waveguiding effect is likely to remain when one considers /c(k, n) rather 
than 7kn- A full answer to this question would require a calculation of V^^ for different k 
and n. 

Since single crystal Fe and Ni already have some intrinsic anisotropy, one might expect 
that this anisotropy could be exploited to obtain a strongly n and k-dependent FOM even 
in single crystals. However, in practice, most magnetic studies of Fe and Ni are carried out 
on polycrystalline samples, which no longer have this anisotropy. The present work provides 
a possible way of recovering this anisotropy, and even more, by use of a periodic lattice of 
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inclusions. 

The present work can be generalized in various other ways. For example, if a homoge- 
neous magnetic layer is perturbed by a periodic array of spin torque oscillators, this would 
generate an artificial magnetic supcrlattice, because the spin torque would provide another 
contribution to He//. Another possibility is to extend the present work to magnonic crys- 
tals with three- dimensional periodicity, though this might be an experimental challenge. The 
present work could conceivably have applications, e. g., in magnonic circuits which exploit 
the strong anisotropy in magnon lifetimes found in the present work. 

In summary, we have calculated the spin wave spectrum of a magnetic supcrlattice with 
two-dimensional periodicity, including for the first time the effects of dissipation. We find a 
striking anisotropy of the spin wave figure of merit, which for typical materials is much larger 
in certain bands near particular points of symmetry in the Brillouin zone. This anisotropy 
implies that propagating spin waves will have much longer lifetimes at certain frequencies 
and in certain directions in k-space , which could be interpreted as a waveguiding effect for 
these excitations. We suggest that this anisotropy might be further increased with suitable 
tuning of the array parameters. 

Funding for this research was provided by the Center for Emergent Materials at the Ohio 
State University, an NSF MRSEC (Award Number DMR-0820414). 
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FIG. 1: (Color online) Left panel: band structure for a triangular lattice of Fe cylinders in Ni, with 
lattice constant a = 10 nm, Fe filling fraction / = 0.5, and no Gilbert damping. Other parameters 
are given in the text. Center panel: same as left panel but for a square lattice. Right panel: 
Brillouin zone for square and triangular lattices with symmetry points indicated. 
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FIG. 2: (Color Online) Left panel: same as center panel of Fig. 1, but with Gilbert damping pa- 
rameters ape = 0.0019 and a^i = 0.064. The widths of the cross-hatched regions are proportional 
to the figm'e of merit (FOM) for the given band, as defined in the text. Right panel: FOM for the 
fourth lowest spin wave band, as a function of superlattice wave vector k, along specified directions 
in the superlattice Brillouin zone (SBZ), and at three diff'erent filling fractions /. 
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FIG. 3: (Color online.) Same as Fig. 2 but for a triangular lattice of Fe cylinders in Ni, with lattice 
constant a = 10 and / = 0.5 (left panel) and / = 0.1, 0.5, and 0.9 (right panel). 
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FIG. 4: (Color Online) Same as Figs. 2 and 3, but showing the FOM as a function of filling fraction 
/ for /ioHo = O.IT, a = 10 nm. 
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